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Abstract 

We investigate a model of a U(l)-Higgs theory on the lattice with compact gauge fields 
but completely suppressed (elementary) monopoles. We study the model at two values of 
the quartic Higgs self-coupling, a strong coupling, A = 3.0, and a weak coupling, A = 0.01. 
We map out the phase diagrams and find that the monopole suppression eliminated the 
confined phase of the standard lattice model at strong gauge coupling. We perform a 
detailed analysis of the static potential and study the mass spectrum in the Coulomb 
and Higgs phases for three values of the gauge coupling. We also probe the existence of a 
scalar bosonium to the extent that our data allow and conclude that further investigations 
are required in the Coulomb phase. 
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1 Introduction 



For many years the compact formulation of QED on a lattice Q has been advocated 
as a preferable alternative to the standard non-compact one. There are at least two 
arguments in favor of quantum electrodynamics with a compact gauge group ||, || . The 
abelian gauge group of QED has to be compact if it appears from some non-abelian, 
unified gauge group. For compact groups the charge is quantized, which is not the case for 
the non-compact theory. Both formulations are supposed to have identical perturbative 
series in the continuum limit, but may lead to completely different infrared behavior on 
a nonperturbative level. 

If we couple U(l) gauge fields to dynamical scalar fields - scalar QED - we obtain 
the simplest prototype of a unified theory. The phase structure of the corresponding 
compact lattice theory and the nature of the vortex string excitations have been intensively 
studied over the years (see, for instance, [^-pT]). Scalar QED shows confinement at strong 
coupling similar to QCD and Coulombic behavior at weak coupling corresponding to the 
real world of electromagnetic phenomena. For sufficiently large values of the parameter 
v - the classical position of the minimum of the Higgs potential - the theory exhibits the 
conventional Higgs behavior. Nevertheless, the nature of the Coulomb-Higgs transition, 
in comparison with the chiral transition, as well as the spectral content of the theory in 
the Coulomb and the Higgs phases deserve further studies. In the non-compact case with 
non-linear Higgs field this has been done recently [O]. The transition turned out to be 
most likely second order and the critical behavior A</> 4 like, i.e., trivial in the continuum 
limit. 

When scalars in the compact theory are in the fundamental representation there is no 
phase boundary that completely separates the confined and Higgs phases. The confined 
and Higgs regions in the phase diagram remain analytically connected. This observation, 
first made for lattice models with frozen out radial modes of the Higgs field ||, was 
later confirmed for more realistic theories with radially varied Higgs fields and for different 
gauge groups. It has led to the conjecture of a principle of complementarity: a confining 
theory of fermions and gauge bosons can be analyzed as if a dynamical Higgs phenomenon 
takes place or as if the gauge symmetry is unbroken [13]. 

Provided all results of lattice calculations can be applied to continuum physics, the 
existence of the Higgs-confinement phase could entail rather strong physical consequences 
(see, e.g., [13, 14 1). However, at finite lattice spacing, lattice theories suffer from lattice 



artifacts. These can be simply scaling violations, i.e., effects from irrelevant operators 
that vanish in the continuum limit, as is believed to be the case for lattice QCD. Or they 
can consist of field configurations on the lattice which have no continuum counterpart 
but strongly influence the infrared properties of the lattice theory. The confinement 
phase in the compact abelian lattice theory is a non-physical phenomenon traced back 
to the particular discretization of gauge fields on the lattice. The confinement-Coulomb 
phase transition is due to the condensation of lattice monopoles, i.e., artifacts [15]. It 
is straightforward to modify Wilson's theory such that monopoles are removed from the 



functional measure in the euclidean functional integral fig , 17 1. 

In this paper we want to discuss the abelian gauge-Higgs theory on the lattice with 
compact U(l) gauge fields in a modified version, with monopoles suppressed. As shown 
m 

U, for the pure U{\) gauge theory the phase transition to the confinement phase 
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disappears at real bare coupling leaving only a unique Coulomb phase [17]. The suppres- 
sion of monopoles in the U(l) theory with staggered fermions entails the disappearance 
of the chiral transition in the zero-mass limit (at least, in the quenched approximation) 
II- 

In the theory with monopoles suppressed the Higgs phase is not anymore part of 
the Higgs-confinement phase, similar to the non-compact formulation, as a confinement 
phase does not exist. The (conjectured) complementarity principle therefore can not be 
applied in this case. It is interesting to see to what extent this influences the physics 
(spectrum, etc.) in the Higgs phase. 

Dealing with the modified theory also has some technical advantages in comparison 
with the standard compact version when considering the Coulomb phase. Since monopoles 
are completely suppressed, problems related to their existence disappear. For instance, it 
is easier to identify the lowest photonic states in plaquette-plaquette correlators, because 
of the smaller overlap with higher energy states with photon quantum numbers [17|. 

Here we discuss the complete phase structure of this model including also imaginary, 
i.e., unphysical bare coupling values. Within the Higgs and the Coulomb phases we try 
to identify and to interpret the lowest lying gauge and Higgs field excitations. This is a 
technically important and, as we shall see, difficult task. 

In the second section we formulate our model and introduce all notations. The third 
section is devoted to the study of the phase structure of the modified theory and its com- 
parison with that of the standard one. In the fourth section we investigate the behavior 
of the heavy charge potential, the screening energy and the Fredenhagen-Marcu order 
parameter in the vicinity of the Higgs phase transition. The fifth section is devoted to 
the study of the spectrum in the Higgs and Coulomb phases. The last section is reserved 
for conclusions. 



2 The Model 



The compact lattice U(l) Higgs model with a scalar field 4> of charge one is defined by 
the action 



SuH 



i 



1--(U D + U*) 



X [1=1 

+A£(0*^-1) 2 + 5>*0 X 



(1) 



where the first term corresponds to the Wilson action Syy with Un as the product of the 
link variables U\ £ U(l) around a plaquette and the remaining terms correspond to the 
A<^ 4 continuum action with minimal coupling to the gauge fields, 



d 4 x[\D^(p x \ 2 + m 2 \<p x \ 2 + \\<Px\% 



(2) 



The relationship between the parameters of the continuum and lattice actions are: 

a<Px = 4>xVk, A = A/k 2 , (am) 2 = (1 — 2A — 8k) /k . (3) 
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4> can be decomposed into a radial and angular part as follows: 



(t>x = PxVx; Px G cr x G U(l). (4) 

The Wilson action SV which is the first term in eq. (|l|) may be given in terms of the 
plaquette angles 9 X -^ U G [— vr, n) 

% = /3^(l-cose x;tlv ). (5) 

X\puV 

where, given the link angles 9 X „ G [— 7r,7r) with £4,,, = e ldx ^ 1 , the plaquette angle is 
expressed as 9 X .^ V = 9 XjfJi + 9 X+ ^ U - 9 x+0 ^ - 6 X>1/ . The physical flux 9 X -^ G [—k, tt) is then 
defined as 

$x\[iv — @x;iiv ~\~ 27r7la;;/iZ/ (6) 

where n x -^ u measures the number of Dirac strings. The modified Wilson action supple- 
mented by a monopole term Sma [01 17] is of the form 

SmA = S w + 7^ \M x . p \ (7) 

x;p 

where 7 is a chemical potential controlling the density of monopoles and M x - p measures 
the net monopole flux out of the 3D cube orthogonal to the direction p: 



2ttM x - p — -e pa ^ u {9 x+ a- P v — 9 x - pv ). (8) 



Using relations @ and (||), eq. (jl]) takes the form 

Sue = P ^2 (I - cos Ox-fa) - Ik PxPx+fi cos(a x+fl - a x + 9 X;/J ) 

xyiv x fi 

+AE(^-!) 2 + E^- (9) 

X X 

The modified compact U(l)-Higgs model, henceforth called modified scalar QED, is de- 
fined by the action Smuh which is just Sjjh with Sw replaced by Sma- 

Smuh = S uh + iY,\ M ^p\- ( 10 ) 

In our simulations, the condition 7 = 00 is set in (|l0|), ensuring that monopoles are 
completely suppressed. Notice that in previous investigations [17, 18] negative plaquette- 
values have been suppressed, too. 



3 Phase Diagram 

We first explored the phase diagram of the model ( |lO| ) with 7 = 00 in the non-linear 
limit A = 00, where the radial part p x of the Higgs field is frozen to unity. This was 
done with "hysteresis" type runs on a 6 4 lattice. We monitored the average plaquette, 
(Re(fTo)), and the 'average link', {(jfUcj)). The phase diagram found in this way is shown 
in Figure |l[ At positive (3 we found the phase transition line separating the Coulomb 
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Figure 1: The phase diagram for A = 0.01 and A = oo. 



phase at small k from the Higgs phase at larger k. Because of the complete suppression 
of monopoles no confining phase was found, and the Higgs transition line did not end - 
in fact could not end - somewhere in the interior of the (/?, At)-plane. This is in contrast 
to the model without suppression of monopoles, where the confinement - Coulomb phase 
transition of the pure gauge model turns into a phase transition line that meets the Higgs 
transition line, and where the extension of the Higgs transition line ends in the interior 
of the (/?, fv)-plane. This ending is possible since before reaching it, the Higgs transition 
lines separated the confinement and Higgs phases which are both massive phases, and are 
analytically connected. As for the unrestricted model, hysteresis runs do not give clear 
evidence for the order of the Higgs transition in the model with completely suppressed 
monopoles, though a weak first order transition seems somewhat more likely. 

Though the model is not physical there since the bare gauge coupling becomes imag- 
inary, one can study the model on the lattice also for negative (3. In our model with 
completely suppressed monopoles, nothing special seems to happen at (3 = 0. The Higgs 
transition continues to negative (3 until (3 ~ —0.7. The pure gauge model has a transition 
at (3* ~ —0.7 into a frustrated phase. The Wilson pure gauge action has a symmetry 
(3 — ► —(3 (with Un — * —Un), but the monopole term prefers positive Un, destroying the 
symmetry and leading to frustration. This first order transition continues into the ((3, k)- 
plane - we have seen long lived coexisting states at ((3,k) = (—0.705,0.15) -, meets the 
Higgs transition and continues upward. Around (f3, k) = (—0.9, 0.65) the transition line 
splits into two, with one branch soon running almost horizontally, parallel to the (3— axis. 
We did not spend any resources to either locate these transitions precisely or study the 
properties of the different phases in the negative (3 half-plane, since we believe that these 
phases are irrelevant for the continuum limit. We note, though, that both the Higgs and 
Coulomb phases appear to continue into the negative j3 region. 

We next explored the phase diagram for a weaker scalar self-coupling, A = 0.01. For 
the model without the suppression of monopoles it is known that the phase diagram 
changes from large to small self-coupling qualitatively Q. While for large A the Higgs 
transition ends in the interior of the (/3, ft)-plane, for smaller couplings it extends to, and 
beyond, the (3 = axis. The phase diagram for our model with complete suppression of 
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monopoles is shown in Figure |l[ It resembles the phase diagram for the non-linear case 
(A = oo): the Higgs transition separating Coulomb and Higgs phase at positive (3, now 
clearly of first order, continues to negative (3 until it is met by the "frustration" transition 
line emanating from (3 ~ —0.7 on the pure gauge axis. Unlike the non-linear case, for 
A = 0.01 two distinct phase transition lines appear to be emerging from this point and to 
continue at different angles to larger k and more negative [3. 

In the remainder of this paper we will concentrate on the physics around the Higgs 
phase transition, and we will present evidence for the Coulomb and Higgs character of 
the phases. Our results are based on between 4000 and 8000 essentially statistically 
independent measurements, after thermalization, for each A, /3, k triplet, except for the 
points at A = 3.0, (3 = —0.5 where we have only 2000 measurements. 



4 Potential, Screening Energy and Order Parameter 

In gauge theories, unlike in spin models or scalar field theories, there do not exist local 
order parameters to distinguish the different phases. Rather one has to study non-local 
objects for this purpose. For pure gauge theories, for example, area or perimeter behavior 
of Wilson loops distinguishes between confining and non-confining phases. However, it 
has been known for some time that the Wilson loop criterion fails in gauge theories where 
dynamical matter fields screen the potential H|. In such theories, gauge-invariant 
objects involving external sources seem to serve better to study the structure of the 
model. 

Therefore, besides computing the usual Wilson loops we also studied the gauge- 
invariant two-point function G(T,R), defined by 

G(T,R) = (<t>* x l[U l( l>y), \x-y\=T (11) 
/er 

where the Ui's are link variables on an oriented path T connecting the points x and y as 
follows: 

V x 



x x- 



R 

Notice the ordering of arguments in G(T, R). The first gives the direction in which one side 
does not contain a string of [/-fields. We also computed the gauge-invariant two-point 
function G(R, T) where the matter fields are in the same time-slice. 

These gauge-invariant functions turn out to be much better objects to study confine- 
ment and/or Higgs character || in compact scalar QED. In our case, they turn out to be 
very good order parameters for signalling the Higgs transition. 

To possibly improve the signal to noise ratio we not only measured Wilson loops 
and the gauge-invariant two-point functions with the ordinary link fields but also with 
"APE-smeared" links [[HJ for the spatial segments. Use of such smeared links has been 
extremely helpful in the study of the heavy quark potential in lattice QCD. 
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4.1 Potential 



To be able to study the static potential we have computed on-axis and some off-axis 
Wilson loops, the latter only for A = 0.01. We considered the distances R = n, and in 
the case of having off-axis loops also y/2n and y5n. To try to get a better overlap with 
the ground-state in 

W(R, T) = cexp{-T • V(R)} + subleading terms (12) 

we computed the Wilson loops also with smeared spatial links. We shall refer to these 
Wilson loops as "smeared" loops. For A = 3.0 and for A = 0.01 at (3 = 2.5, actually only 
smeared loops were computed. In simulations of lattice QCD the use of smeared loops 
leads to an early plateau in the finite T approximants (effective potentials) 

V T (R) = loJ W ^ T) ) (13) 

to the potential V(R) and is invaluable in the extraction of the potential. 

For unsmeared planar Wilson loops, using the fact that they are symmetric under the 
exchange of space and time direction, we can estimate how fast the finite T approximants 
approach the true potential. For a Yukawa potential we find 



V T (R) = V(R) + cR 



ie^ T - 1 e -M(T+i) 
T T+l 



(14) 



i.e., an exponential approach. For a Coulomb potential, on the other hand, the \i — ► 
limit of (|14|), the approach is only power-like, like 1/T 2 . 

For the U(l)-Higgs model studied by us, smeared Wilson loops were of some use, 
though not nearly as dramatic as in QCD. They lead to a slightly earlier plateau in the 
Higgs phase, and to smaller errors in both Higgs and Coulomb phases. In the Coulomb 
phase, the smeared loops gave a somewhat smaller dependence of Vt{R) on T, but they 
still seemed to approach an asymptotic value only as 1/T 2 , just like the potential approx- 
imants obtained from normal Wilson loops. We generally did get a good signal even at 
the largest distances R, out to T = 7, on the 8 3 x 16 lattices that we simulated. 

For pure gauge U(l) the perturbative expansion of the lattice potential leads to 



V^\R)=g 2 v£ on \R) 



g 2 

1 + T 



+ 0(g 6 ), (15) 



where the finite T approximants are defined analogous to (|i~3|), g 2 = 1/(5 and V>p oul (R) is 
the tree-level - T dependent - contribution with the coupling factored out. Possible non- 
Coulomb terms appear only at higher orders, and to order g 4 we just find a renormalization 
of the bare coupling. 

Scalar fields, as long as we stay in the Coulomb phase of course, induce some non- 
Coulomb terms already at order g 4 . They are, however, except possibly very close to 
the Higgs transition, expected to be small. Since, as mentioned before, the finite T 
approximants to the potential approach their asymptotic value, the true potential, only 
as 1/T 2 , we decided, rather than trying to extract this 'true' potential and then attempting 
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Table 1: The couplings g R in the Coulomb phase from fits as in eq. (|i~6|) of the finite T 
approximants to the potential to the leading finite volume lattice perturbative values for 
A = 0.01. Notice that the perturbative value for (3 = 2.5 is 0.44 plus O(0. 4 3 ). 

fits to it, to make fits of the finite T approximants to the lowest order perturbative from, 
but with a renormalized coupling as the (only) fit parameter, 

V^ c (R)=g 2 R V^\R) . (16) 

To fit also the effective potential from smeared Wilson loops, we computed the 'smeared' 
analogue of Vp onl (R) with the smeared links, constructed as in the Monte Carlo simu- 
lations, but expanded to lowest order in perturbation theory. These fits work very well 
(for positive (3) once we exclude the approximants from T = 1 and the shortest distance, 
R = 1. The results are listed in Table [l] for A = 0.01 and in Table ^ for A = 3.0. To 
illustrate the quality of the fits we show in Figure || the difference between the measured 
and fitted effective potentials for T = 1, which was not included in the fit, and for T = 4 
and 7, for (3 = 0.5, A = 0.01 and k = 0.1475. For j3 = 2.5, the bare coupling is g 2 = 0.4 
and we expect perturbation theory to be applicable. Indeed, the difference between the 
values obtained from the fits in Tables |l| and § to the perturbative prediction of 0.44 from 
(|l5| ) is certainly compatible with coming from higher order perturbative corrections. 

As can be seen from Table [l| the results from normal and smeared loops agree within 
errors with the smeared loops giving somewhat smaller errors. The dependence on k, i.e., 
the scalar fields, is rather weak with g 2 R decreasing slightly as the scalars become lighter 
(as the Higgs transition is approached) as one would expect. 

In Table [l] no g R is quoted for normal Wilson loops for (3 = —0.5. The reason is 
that several loops, including the plaquette, were negative there. In fact, exactly those 
of the small loops were negative, that would become positive under the transformation 
that makes the symmetry (3 — ► —/3 for the Wilson model, i.e., without the suppression 
of monopoles. Smearing got rid of this problem. However, the Wilson loops remained 
rather noisy compared to those obtained at positive gauge coupling and the approach to 
an asymptotic value of the effective potentials ( |l~3|) was not monotonia This fact points to 
a non-positive transfer matrix, as one might expect for the imaginary bare gauge coupling 
one has for negative (3. A similar, non-monotonic behavior was observed in the Higgs 
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Table 2: The couplings g\ in the Coulomb phase from fits as in eq. (^) of the finite T 
approximants to the potential to the leading finite volume lattice perturbative values for 
A = 3.0. Notice that the perturbative value for (3 = 2.5 is 0.44 plus 0(0. 4 3 ). 
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Figure 2: The difference between the measured and fitted effective potentials for T = 1, 
4 and 7, for (3 = 0.5, A = 0.01 and k = 0.1475. The values for T = 4 and 7 have been 
displaced by 5t = 0.1 and 0.2 vertically for better viewing. 
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Table 3: The fit parameters from fits to the finite T approximants to the lattice Yukawa 
potential, obtained from lowest order Wilson loops with a massive gauge boson propaga- 
tor, for A = 0.01. The first line for each data set is the fit with smallest x 2 / n df, and the 
second line comes from T = 7, unless the fits are identical. 



phase, as will be discussed further below. We should add here, that contrary to all other 
cases, only time separations up to T = 4 could be included in the fit for g^. For larger T 
we did not get a signal for the effective potential except for the smallest distances. 

In the Higgs phase, the potential is a screened, or Yukawa potential. There we made 
fits to finite T approximants of the lattice Yukawa potential 

V™ c (R)=g 2 R -V™(R) , (17) 

obtained from lowest order Wilson loops with the vector boson propagator in gauge 
with £ = 1 given by 

G,iu{p) = j=; — . 2 , a ■ u2/h' ( 18 ) 
4 2^ sin + 4 smh h£ 

with fiy the mass of the gauge boson. We used as a fitting parameter m 7 = 2sinh4f . 

These fits, tabulated in Table [| and ||, worked well for the A = 0.01, j3 = 0.5 data where 
niy is relatively large and a plateau is achieved quite soon, but they worked much better 
for the data at the larger (3 value, (3 = 2.5. A sample fit is shown in Figure ||. Notice the 
violations of rotational symmetry, especially at short distances, which are well followed by 
the lattice Yukawa potential. A fit to the continuum Yukawa potential -^e -7 ™ 7 ^ would not 
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Table 4: The fit parameters from fits to the finite T approximants to the lattice Yukawa 
potential, obtained from lowest order Wilson loops with a massive gauge boson propaga- 
tor, for A = 3.0. 



have worked. For data in the Coulomb phase, the Yukawa fits, correctly gave a vanishing 
photon mass, within error. 

It maybe worthwhile to note that the attempt to fit to the 'naive' lattice Yukawa 
potential, the T -> oo limit of V^ uk (R), 

VnaUR) =V + &± £ 4V ^^Stt 2 ( 19 ) 

was successful only for A = 0.01 and /3 = 0.5. For (3 = 2.5 and both values of A, 
as well as for f3 = 0.5 and A = 3.0 we could not obtain fits with a decent x 2 using 



eq. (19). This observation is in disagreement with some statements in the literature ||. 
A possible explanation might be that smaller statistical errors in our case lead to a better 
discrimination between different fit formulae, thus excluding the 'naive' fits. 
In the quasiclassical perturbative approach one expects 

m 2 - — - <0V>- (20) 



Perturbative arguments suggest that this equality should be valid in the weak coupling 
region when the dimensionless mass /U 7 is sufficiently small, It is 

interesting to note that this equality seems to hold even in the region where m 7 is not 
very small. In Figure || we compare the values of m 7 taken from Table || with that 



calculated via eq. (20) at (3 = 0.5 and A = 0.01. Qualitatively, the agreement is quite 
reasonable, indicating consistency of the quasiclassical estimation. 

For (3 = —0.5 the potentials looked rather strange, with the point at R = 1 being the 
highest. An example is shown in Figure |5|. The potentials from smeared and normal Wil- 
son loops, where available, agreed. As in the Coulomb phase, the potential approximants 
show a non-monotonic behavior with T, though a plateau is reached rather quickly. Even 
ignoring the R = 1 point, no good fits could be obtained, and hence no conclusions drawn. 
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Figure 3: Monte Carlo potential from smeared Wilson loops and the fit to the effective po- 
tential extracted from lowest order perturbative Wilson loops with a massive propagator, 
for = 2.5, A = 0.01 and k = 0.1375. 
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Figure 4: The parameter m 7 in the Higgs phase. Diamonds correspond to m 7 values 
taken from Table ||. The solid line corresponds to r.h.s. of eq. (|20|). 
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Figure 5: Monte Carlo potential from smeared Wilson loops for j3 
K = 0.1875. 
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4.2 Screening Energy 

The gauge invariant two-point function ( pd| ) can be interpreted as an external static 
charge propagating forward in (Euclidean) time. The presence of the external charge 
manifests itself as an energy increase fx of the lowest state of the system. The screening 
of an external source may be studied from the exponential fall-off of G(T, R) for fixed R 
and for T large enough: 

G(T, R) ~ f(R) exp(-/xT). (21) 

For smaller T there are subleading exponentials to this, /i is the lowest energy in the 
fields that screen an external charge, and is therefore called the screening energy. 

In a confining phase in the presence of matter fields or in the Higgs phase, the external 
charges are always completely screened and one therefore finds that 

H = iy(oo) (22) 

with V{R) the static potential that can be extracted from Wilson loops. In the Coulomb 
phase the screening energy is the energy of a bound-state of a light charge, represented by 
the field 4>* in ( |TT|) and the static charge, corresponding to the straight time-like segment. 
With m c the mass of the light charged particle the screening energy is then 

H = -y(oo) + m c - E b (23) 

with Eb the binding energy. One therefore expects [i > V(oo)/2 in the Coulomb phase. 

The gauge-invariant two-point function G(T, R) could be calculated quite accurately 
and showed in both the free charge phase and the Higgs phase the expected exponential 
falloff as a function of T, for T large enough with R fixed as in eq. (21). Figure || shows 



the behavior of the two-point function as a function of T and different R for k = 0.1506, 
(3 = 0.5 and A = 0.01 near the phase transition. The screening energy, as extracted 
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P, n\fj, 


R = 1 


R = 2 


R = 3 


R = 4 


±V(R = 4) 


-0.5, 0.1850 


0.1635(1) 


0.1634(1) 


0.1632(3) 


0.1630(3) 


0.1637(10) 


-0.5, 0.1875 


0.1325(1) 


0.1327(1) 


0.1330(2) 


0.1329(2) 


0.1322(5) 


-0.5, 0.1900 


0.1123(1) 


0.1128(2) 


0.1128(2) 


0.1129(3) 


0.1117(3) 


0.5, 0.1450 


1.25(16) 


0.94(25) 


1.37(32) 


1.60(10) 


0.34(1) 


0.5, 0.1475 


1.00(5) 


0.76(7) 


0.68(35) 


0.57(49) 


0.32(1) 


0.5, 0.1500 


0.63(5) 


0.65(5) 


0.59(10) 


0.47(10) 


0.33(1) 


0.5, 0.1506 


0.2437(5) 


0.2442(6) 


0.2447(8) 


0.2412(10) 


0.242(3) 


0.5, 0.1525 


0.1844(2) 


0.1841(2) 


0.1842(3) 


0.1842(4) 


0.183(1) 


0.5, 0.1550 


0.1484(1) 


0.1486(1) 


0.1487(2) 


0.1485(3) 


0.1487(6) 


0.5, 0.1575 


0.1264(1) 


0.1265(1) 


0.1264(1) 


0.1263(1) 


0.1267(3) 



Table 5: Screening energy n for A = 0.01 extracted from fits with unsmeared operators 
for different R values and k, generally with T > 4. 



P, k\h 


R = l 


R = 2 


R = 3 


R = A 


~ 2 V(R = 4) 


-0.5, 0.1850 


0.1635(1) 


0.1637(1) 


0.1637(1) 


0.1637(1) 


0.1634(2) 


-0.5, 0.1875 


0.1324(1) 


0.1325(1) 


0.1326(1) 


0.1326(1) 


0.1323(1) 


-0.5, 0.1900 


0.1124(1) 


0.1124(1) 


0.1122(1) 


0.1124(1) 


0.1123(1) 


0.5, 0.1450 


1.31(13) 


0.73(12) 


1.49(7) 


1.37(7) 


0.332(1) 


0.5, 0.1475 


1.05(6) 


0.98(6) 


1.00(9) 


1.05(10) 


0.332(2) 


0.5, 0.1500 


0.64(5) 


0.65(5) 


0.49(4) 


0.64(7) 


0.329(2) 


0.5, 0.1506 


0.2435(5) 


0.2438(6) 


0.2438(6) 


0.2430(6) 


0.2426(7) 


0.5, 0.1525 


0.1841(2) 


0.1840(3) 


0.1842(3) 


0.1845(2) 


0.1838(3) 


0.5, 0.1550 


0.1484(1) 


0.1484(1) 


0.1485(1) 


0.1485(2) 


0.1485(2) 


0.5, 0.1575 


0.1263(1) 


0.1264(1) 


0.1264(1) 


0.1264(1) 


0.1263(2) 



Table 6: Screening energy fi for A = 0.01 extracted from fits with smeared operators for 
different R values and k, with T > 4. 



14 



2 



R = 1 H^H 
R = 2 ^ 
R = 3 h^h 
R = 4 £X— l 

m 



1.5 



a* 

H 

s 



0.5 



-0.5 




Figure 6: The logarithm of G(T, i?) plotted as a function of T for various R at k = 0.1506, 
/? = 0.5 and A = 0.01 . 



from different distances R, is tabulated in Tables || and |6] for unsmeared and smeared 
operators, respectively. The last column in both Tables is an estimation of the selfenergy 
of an isolated charge E q = \V(R = oo) from the potential V(R = 4, T = 5). The 
screening energy as defined by eq. ([HI) is supposed to be independent of the distance R in 
the gauge invariant two-point function G(T, R) B. Our results are in general agreement 
with this expectation. In the Coulomb phase, where /U is rather large and hence less 
well determined since G(T, R) becomes noisier at larger T, there are some variations, but 
clearly no definite trend. We also note that both smeared and unsmeared operators work 
equally well leading to consistent values of fj,, except for the 'problem points' mentioned 
above. We therefore expect that with much better statistics [i would again turn out to 
be independent of R. 

In Figure [5], we show the screening energy [i as a function of k, extracted from fits 
to log(G(T, R)). Also shown in the figure is V(R = 4)/2, our approximation to E q = 
V(oo)/2, the selfenergy of an isolated external charge. In the Higgs phase, V(R = 4)/2 
is expected to be an excellent approximation to V(oo)/2, and as expected \i and V(oo)/2 
agree very well, (^). In the Coulomb phase, V(oo)/2 is slightly higher than V(R = 4)/2 
but, as expected, the screening energy fj, is clearly much larger. 

For A = 3.0 and all /3-values, and for A = 0.01 at (3 = 2.5, due to an oversight we only 
computed G(R,T) used for the Fredenhagen-Marcu order parameter, but not G(T,R). 
Therefore we have no data for the screening energy for those parameter values. 



15 



1.6 



1.4 



1.2 



II. s 



0.6 



0.4 



0.2 



screening energy 
selfenergy 



0.144 



0.146 



0.148 



0.15 0.152 
Kappa 



0.154 



0.156 



0.158 



0.16 



Figure 7: The screening energy p and the selfenergy of an isolated external charge, E q , 
approximated by V(R = 4)/2, as a function of k for (3 = 0.5 and A = 0.01. p is determined 
from the exponential decay of G(T, R = 1). 



4.3 Fredenhagen-Marcu Order Parameter 

Fredenhagen and Marcu found that a suitable criterion to study confinement in gauge 
theories with matter fields |21] can be formulated in terms of the ratio 



P fm(R,T) 



This ratio has the asymptotic behavior 



Pfm(R,T) 



R — *oo,T—R/2 







G(R,T) 
W{R, 2T)V2' 



confined or Higgs phase 
Coulomb phase 



(24) 



(25) 



The physical interpretation is that one of two charges is removed as R — » oo. If free 
charges can exist the overlap of the remaining one-charge state with the vacuum is zero 
and the ratio (24) vanishes. If the remaining charge is screened as in a confined or Higgs 
phase the ratio attains a finite, non-zero limit. In both cases the division by the Wilson 
loop ensures the cancellation of the perimeter law behavior in G(R,T). The asymptotic 
value Pf'm is therefore called the vacuum overlap order parameter. The ratio between T 
and R is held fixed in the limiting process to ensure projection onto the state with lowest 
energy. 

The Fredenhagen-Marcu parameter pfm{R, T) turns out to be a very good parameter 
signalling the Higgs phase transition H. As our lattices were restricted to 8 3 x 16, we 
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Figure 8: The Fredenhagen-Marcu parameter pfm(R,T) as a function of k for /? = 0.5 
and A = 0.01. 



a? 




Figure 9: The Fredenhagen-Marcu parameter pfm(R,T) as a function of k for (5 = 2.5 
and A = 0.01. 
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Figure 10: The Fredenhagen-Marcu parameter pfm(R,T) as a function of k for f3 = 0.5 
and A = 3.0. 
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Figure 11: The Fredenhagen-Marcu parameter pfm{R,T) as a function of k for (3 = 2.5 
and A = 3.0. 
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P, K\ PFM 


R = 2, norm 


R = 4, norm 


R = 2, smear 


R = 4, smear 


-0.5, 0.1850 


0.278(2) 


0.007(6) 


0.0542(4) 


0.0030(8) 


-0.5, 0.1875 


21.137(5) 


13.886(5) 


15.082(3) 


13.015(3) 


-0.5, 0.1900 


19.12(4) 


19.19(6) 


17.308(3) 


15.279(3) 


0.5, 0.1450 


0.0654(2) 


0.0024(3) 


0.0420(1) 


0.0025(1) 


0.5, 0.1475 


0.0786(3) 


0.0041(3) 


0.0510(2) 


0.0045(2) 


0.5, 0.1500 


0.091(4) 


0.020(3) 


0.102(5) 


0.036(6) 


0.5, 0.1506 


1.89(2) 


1.93(2) 


1.88(2) 


1.91(2) 


0.5, 0.1525 


4.47(1) 


2.673(9) 


4.148(9) 


4.19(1) 


0.5, 0.1550 


5.556(5) 


4.652(5) 


6.217(5) 


6.258(6) 


U.D, U.lCWO 


a arte. ( n\ 
O.OOO^Z J 


O.Z04(Z J 


I .\)lo\6 J 


1 .\)OZ{6 ) 


2.5, 0.1275 






0.0433(2) 


0.0034(3) 


2.5, 0.1300 






0.0530(4) 


0.0080(6) 


2.5, 0.1325 






0.36(1) 


0.30(1) 


2.5, 0.1350 






1.536(4) 


1.424(3) 


2.5, 0.1375 






2.654(2) 


2.496(2) 



Table 7: The Fredenhagen-Marcu parameter pfm{R,T) for A = 0.01 extracted from 
unsmeared (normal) and smeared operators for R = 2 and R = 4 with T = R/2. 



were able to calculate this only up to pfm(4, 2). But as Figures |§] to 11 indicate, this 
quantity behaves like an order parameter. 

Our data for A = 0.01 are tabulated in Table ^ from computations with normal and 
smeared space-like links and in Table |8| for A = 3.0 with smeared space-like links. 



5 Spectrum Calculations 

To probe the spectrum of the model, we constructed correlation functions of the following 
operators: 

S (x,t) = ReJ2€U x ,^ x+ft (26) 

p- 

Ov(x,t) = Tm^lU^+n (27) 

P (x,t) = ImJ2 u a (28) 

□ 

with the subscripts S, V, and P corresponding to the scalar boson, the vector boson 
and the photon, and /i, v = 1,2,3. These states have J PC quantum numbers ++ , 1 
and 1^ respectively. The operator Op at momentum p = has the wrong parity for a 
photon state. However, at non-zero momenta, p > 0, parity is no longer a good quantum 
number for extended objects like the plaquettes, Un, and Op will have non-zero overlap 



with photon states [22] 



To study the spectrum of the model, we extract the connected correlations of a generic 
operator O with quantum numbers J PC and momentum p: 
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P, k\p F m 


R = 2 


R = 4 


-0.5, 0.3150 


0.122(2) 


0.057(7) 


-0.5, 0.3175 


0.135(4) 


0.06(1) 


-U.o, U.ozUU 


U.LI L[0 ) 


U.llD^O ) 


-0.5, 0.3225 


0.235(5) 


0.28(1) 


-0.5, 0.3250 


0.287(3) 


0.311(7) 


-0.5, 0.3275 


0.335(3) 


0.371(6) 


0.5, 0.2100 


0.0574(3) 


0.0132(6) 


0.5, 0.2125 


0.0665(8) 


0.021(1) 


0.5, 0.2150 


0.090(2) 


0.049(2) 


0.5, 0.2175 


0.1465(8) 


0.118(1) 


0.5, 0.2200 


0.1913(7) 


0.171(1) 


0.5, 0.2225 


0.2266(5) 


0.2126(8) 


0.5, 0.2250 


0.2610(4) 


0.2504(5) 


2.5, 0.1750 


0.0463(5) 


0.0090(9) 


2.5, 0.1775 


0.058(1) 


0.020(2) 


2.5, 0.1800 


0.103(2) 


0.071(2) 


2.5, 0.1825 


0.1465(7) 


0.1186(9) 


2.5, 0.1850 


0.1823(4) 


0.1574(5) 


2.5, 0.1875 


0.2149(3) 


0.1921(3) 



Table 8: The Fredenhagen-Marcu parameter pfm(R,T) for A = 3.0 extracted from 
smeared operators for R = 2 and R = 4 with T = R/2. 
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C(r) = (0(t + r;p)OHt;p}) - (0(t + r;p))(OHt;p)); 

(29) 

0(t;p) = Y,e^-0{t;x). 

X 

These correlations are a sum of exponentials exp(— EiT) with E{ the energies of the 
different states with the given quantum number. We are interested in the state with 
the lowest energy, Eq, in each channel. This lowest energy can be approximated by the 
effective energy 

Eeff = ~ ln W~^y (30) 

at large enough r with r <C Lt/2. In practice, the situation is more complicated as the 
signals for larger time-slices die away rapidly. The true energy E could be quite different 
from the calculated energy E e f t if there is a lot of contamination from the higher energy 
states. In order to maximize the overlap of the lowest lying state with our operator and 
obtain a better signal for the effective energies, we used the APE smearing technique and 
applied between five and ten smearing iterations. In our computations, we found that we 
were able to get reliable signals up to r = 2. Once the energies E have been determined, 
the masses of the corresponding particles may be calculated through the lattice dispersion 
relation: 



sinh 2 - — = sinh 2 h sin 2 (31) 

2 2^2 



The momenta, here, are given by pi = 2itki/(aL), with k{ = —\L + 1, . . . , ^L. We 
shall refer to these momenta simply as a pi = hi" or as "p = k" , keeping in mind that the 
momenta are given in units of 2ir/(aL). 

5.1 Photon 

The operator Op, which at non-zero momentum couples to the photon, was measured at 
several non-zero momenta. The results for the effective energy at different (3 and A are 
tabulated in Tables ^ to 11. As an example, we show in Figure ^ the effective energies 



for 3 different momenta for A = 0.01 and f3 = 0.5. The lines correspond to the effective 
energies that would be expected for a massless particle at these 3 momenta. As expected, 
the photon is massless at k's smaller than k c (in the Coulomb phase) and acquires a mass 
at k > k c (in the Higgs phase). This mass, extracted from the energy with the dispersion 
relation (j3l|), is given in the last three columns of Tables ^ to 11. We note that the masses 



extracted from the dispersion relation are in general agreement with, though somewhat 
smaller than, those from the Yukawa fits to the potential (see Tables |3| and ||), but they 
tend to have larger errors. The dependence of the photon mass on k above the phase 
transition (PT) is rather strong for small self-coupling, i.e., A = 0.01, and a discontinuity 
of the mass can not be excluded. On the contrary, for large self-coupling (A = 3.0) the 
photon mass grows very slowly above the PT, in agreement with a previous observation 
in |3 . 
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Figure 12: Effective Energies of the 1 H state at three different non-zero momenta as 
indicated for A = 0.01 and (3 = 0.5. The lines correspond to the effective energies expected 
from a massless particle. 

The study of other order parameters in the previous sections convinced us that a 
massless state should exist also at negative /9's (/? > —0.7). This is indeed the case. As an 
illustration we show in Figure [l^ the effective energies for different momenta for A = 3.0 
and = -0.5. 



5.2 Scalar Boson 

It is rather plausible that in the Coulomb phase the scalars <j) interacting through the 
attractive Coulombic force and the </> 4 coupling can form massive bound-states, bosonium, 
analogous to positronium in normal (fermionic) QED. One of the possible bound-states 
has the same quantum numbers, ++ , as the Higgs boson in the Higgs phase. The analogy 
with positronium makes it interesting to subject this state to a somewhat detailed study. 
Therefore, we attempted to check its lattice dispersion relation at different (3 and A, 
including the negative value j3 = —0.5. 

We measured at two different momenta, p = and p = 1. We were especially interested 
to see whether a two free particles dispersion (DR) relation was to be favored over a 
single-particle dispersion, indicating a bound-state, in the Coulomb phase. To investigate 
this, we calculated the mass of the two assumed free particles from the p = energy 
(E p=0 = 2m) and used this mass to compute the energy Ef~ a of the combined two- 
particle p = 1 state, 



Etwo = m + (2/a)arsinhw sinh 2 — — h sin 2 — . 



(32) 
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Figure 13: Effective Energies of the l" 1 state at three different momenta indicated for 
A = 3.0 and (3 = —0.5. The solid lines correspond to the effective energies expected from 
a massless particle. 



(3,K 


Ek=i 


Ek=2 


Ek=3 


m k =i 


m k =2 


m k =3 


-0.5, 0.1850 


0.74(2) 


1.33(5) 


1.78(8) 








0.5, 0.1450 


0.75(1) 


1.29(4) 


1.57(8) 








0.5, 0.1475 


0.75(1) 


1.35(4) 


1.72(7) 








0.5, 0.1500 


0.756(6) 


1.30(5) 


1.72(11) 








0.5, 0.1506 


1.15(3) 


1.51(5) 


2.00(13) 


0.91(4) 


0.84(11) 


1.3(3) 


0.5, 0.1525 


1.64(11) 


2.22(17) 


2.75(49) 


1.51(13) 


1.97(22) 


2.5(6) 


0.5, 0.1550 


2.01(13) 


1.80(14) 


1.80(20) 


1.93(14) 


1.38(22) 


0.9(6) 


0.5, 0.1575 


2.14(14) 


1.94(15) 


3.09(97) 


2.07(15) 


1.59(22) 


2.9(1.2) 



Table 9: The effective energies and, in the Higgs phase, the extracted masses of the photon 
from unsmeared operators at time-slice 2 for A = 0.01 and for three different values of 
the momentum: k = (0, 1,0), k = (0,2,0) and k = (0,3,0). 
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/3,K 


Ek=i 


Ek=2 


t E k =3 


m k =i 


m k =2 


m k =3 


-0.5, 0.1850 


0.72(1) 


1.32(9) 










-0.5, 0.1875 


0.68(11) 


1.30(49) 










0.5, 0.1450 


0.749(8) 


1.28(3) 


1.66(6) 








0.5, 0.1475 


0.75(1) 


1.34(3) 


1.71(9) 








0.5, 0.1500 


0.754(7) 


1.30(3) 


1.66(8) 








0.5, 0.1506 


1.18(2) 


1.55(4) 


1.90(13) 


0.95(3) 


0.93(8) 


1.1(3) 


0.5, 0.1525 


1.62(4) 


1.95(9) 


2.19(19) 


1.49(5) 


1.60(13) 


1.7(3) 


0.5, 0.1550 


2.03(10) 


1.82(9) 


1.90(14) 


1.95(11) 


1.41(14) 


1.1(3) 


0.5, 0.1575 


2.06(11) 


1.95(10) 


2.85(33) 


1.98(12) 


1.60(14) 


2.6(4) 


2.5, 0.1275 


0.758(8) 


1.33(4) 


1.60(8) 








2.5, 0.1300 


0.75(1) 


1.33(4) 


1.44(5) 
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0.6(4) 
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0.84(2) 


1.30(4) 


1.63(7) 


0.40(5) 






2.5, 0.1375 


0.90(2) 


1.40(5) 


1.78(7) 


0.52(4) 


0.54(17) 


0.8(2) 



Table 10: The effective energies and, in the Higgs phase, the extracted masses of the 
photon from smeared operators at time-slice 2 for A = 0.01 and for three different values 
of the momentum: k = (0, 1, 0), k = (0, 2, 0) and k = (0, 3, 0). 
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0.75(1) 


1.22(2) 


1.85(8) 


0.06(14) 




1.0(2) 


2.5, 0.1850 


0.75(1) 


1.27(4) 


1.55(6) 


0.06(14) 






2.5, 0.1875 


0.77(1) 


1.43(5) 


1.56(7) 


0.19(4) 


0.64(15) 





Table 11: The effective energies and, in the Higgs phase, the extracted masses of the 
photon from smeared operators at time-slice 2 for A = 3.0 and for three different values 
of the momentum: k = (0, 1, 0), k = (0, 2, 0) and k = (0, 3, 0). 
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Energy Measured 
Energy Calculated 



0.174 0.176 0.178 0.18 0.182 0.184 0.186 0.188 0.19 

Kappa 

Figure 14: Measured and calculated Energies of the ++ state at p = 1 assuming a 
two-particle dispersion relation for A = 3.0 and (3 = 2.5. 



We then compared the measured energy E 1 ,^ 1 with the calculated energy Ef wo as well 
as with the calculated energy E p s ~^ g based on the one-particle dispersion relation. The 



results are tabulated in Tables 12 to 14, and shown in Figures 14 to [19]. 

In the Higgs phase the two-particle dispersion relation definitely fails, while the single- 
particle dispersion relation is quite consistent for the two different momenta p we have 
chosen. This is, of course, consistent with the notion that the Higgs boson, in our model, 
is an elementary particle. 

However, in the Coulomb phase the situation is less clear. Neither the single nor 
the two-particle dispersion relations work satisfactorily in the case of the strong self- 
interaction (A = 3.0). As an example, we show in Figures 14 and ^ the measured and 
calculated energies with two-particle DR and one-particle DR, respectively. The failure 
of the one-particle DR becomes even more pronounced for negative /3's (see Figures 
16 and |l7| ). A possible explanation could be that, in the Coulomb phase, our energy 
measurements are more contaminated by higher excitations. Unfortunately our signal 
does not allow us to go to larger r to settle this issue. 

In the case of weak self-interaction the one-particle DR looks somewhat preferable, 



as compared to the two-particle DR. In Figures 18 and 19 we show the measured and 
calculated energies with two-particle DR and one-particle DR, respectively, for /3 = 0.5 
and A = 0.01. A qualitatively similar picture was observed for other /3-values for the 
same A. 
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0.184 0.186 
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Figure 15: Measured and calculated Energies of the ++ state at p = 1 assuming a 
single-particle dispersion relation for A = 3.0 and (3 = 2.5. 



Energy Measured 
Energy Calculated 




Figure 16: Measured and calculated Energies of the ++ state at p = 1 assuming a 
two-particle dispersion relation for A = 3.0 and (3 = —0.5. 
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Energy Measured 
Energy Calculated 



0.314 0.316 0.318 0.32 0.322 0.324 0.326 0.328 0.33 

Kappa 



Figure 17: Measured and calculated Energies of the ++ state at p = 1 assuming a 
single-particle dispersion relation for A = 3.0 and (3 = —0.5. 



Energy Measured 
Energy Calculated 



0.144 0.146 0.148 0.15 0.152 0.154 0.156 0.158 0.16 

Kappa 



Figure 18: Measured and calculated Energies of the ++ state at p = 1 assuming a 
two-particle dispersion relation for A = 0.01 and (3 = 0.5. 
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Energy Measured 
Energy Calculated 
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0.144 0.146 0.148 0.15 0.152 0.154 0.156 0.158 0.16 

Kappa 

Figure 19: Measured and calculated Energies of the ++ state at p = 1 assuming a 
single-particle dispersion relation for A = 0.01 and (3 = 0.5. 



0, K \ E e ff 


E p=0 


E p=1 

^meas 


^two 


E p=i 

sing 


-0.5, 0.1850 


1.65(11) 


1.87(14) 


1.91(9) 


1.76(10) 


-0.5, 0.1850(*) 


0.93(3) 


1.23(3) 


1.33(2) 


1.16(2) 


-0.5, 0.1875 


1.06(4) 


1.35(4) 


1.43(3) 


1.26(3) 


-0.5, 0.1900 


1.17(12) 


1.48(43) 


1.52(10) 


1.35(10) 


0.5, 0.1450 


1.84(17) 


2.01(16) 


2.07(15) 


1.93(15) 


0.5, 0.1475 


1.64(10) 
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0.5, 0.1500 


0.42(9) 


1.17(5) 


0.99(6) 


0.85(4) 


0.5, 0.1506 


0.33(2) 


0.81(2) 


0.93(1) 


0.813(7) 


0.5, 0.1525 


0.66(3) 


0.97(2) 


1.14(2) 


0.98(2) 


0.5, 0.1550 


0.91(3) 


1.15(3) 


1.32(2) 


1.15(2) 


0.5, 0.1575 


1.07(5) 


1.25(4) 


1.44(4) 


1.27(4) 



Table 12: The effective energies of the ++ scalar boson at momentum p = and p = 1 
from unsmeared operators at time-slice 2 for A = 0.01. Tabulated are the effective energies 
which were actually measured and for p = 1 also those computed from a two-particle and 
from a single-particle dispersion relation for the given E p=0 . The (*) data corresponds 
to a run from a new start, starting from k = 0.1875. 
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Table 13: The effective energies of the ++ scalar boson at momentum p = and p = 1 
from smeared operators at time-slice 2 for A = 0.01. Tabulated are the effective energies 
which were actually measured and for p = 1 also those computed from a two-particle and 
from a single-particle dispersion relation for the given E p=0 . The (*) data corresponds 
to a run from a new start, starting from n = 0.1875. 



29 





J?P=0 

^meas 


F,P =l 

^meas 


e p=l 

^two 


sing 


-0.5, 0.3150 


0.80(12) 


1.83(17) 


1.24(9) 


1.07(8) 


-0.5, 0.3175 


0.60(8) 


1 /IO/l o\ 

1.43(12) 


1.10(5) 


n r\ a ( cr \ 

0.94(5) 


-0.5, 0.3200 


0.52(4) 


1.23(9) 


1.05(3) 


0.90(2) 


-D 5 D ^995 


D 61 (*\\ 


1 05(71 


1 1 1 C\\ 

± . ± J- t KJ I 


D 

VJ. cfO I O ) 


-0.5, 0.3250 


0.76(5) 


1.14(5) 


1.21(3) 


1.04(3) 


-0.5, 0.3275 


0.74(5) 


1.04(4) 


1.19(3) 


1.03(3) 


0.5, 0.2100 


1.24(11) 


1.43(9) 


1.58(9) 


1.41(9) 


0.5, 0.2125 


0.71(4) 


1.30(8) 


1.18(3) 


1.01(2) 


0.5, 0.2150 


0.47(2) 


1.16(4) 


1.02(1) 


0.87(1) 


0.5, 0.2175 


0.58(2) 


1.15(5) 


1.08(1) 


0.93(1) 


0.5, 0.2200 


0.74(3) 


1.14(3) 


1.20(2) 


1.03(2) 


0.5, 0.2225 


0.81(3) 


1.17(4) 


1.25(2) 


1.08(2) 


0.5, 0.2250 


0.94(3) 


1.20(4) 


1.34(2) 


1.17(2) 


2.5, 0.1750 


1.07(6) 


1.54(8) 


1.44(5) 


1.27(5) 


2.5, 0.1775 


0.58(5) 


1.32(10) 


1.08(3) 


0.93(3) 


2.5, 0.1800 


0.49(4) 


1.07(4) 


1.03(2) 
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Table 14: The effective energies of the ++ scalar boson at momentum p = and p = 1 
from smeared operators at time-slice 2 for A = 3.0. Tabulated are the effective energies 
which were actually measured and for p = 1 also those computed from a two-particle and 
from a single-particle dispersion relation for the given E p=0 . 
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Figure 20: Effective Energies of the 1 vector boson at p = 1 and p = for A = 0.01 
and = 0.5. 



5.3 Vector Boson 

Correlations of the operator Oy allow, in principle, the measurement of the 1 vector 
boson mass. This works reasonably well in the Higgs phase. As an example, we show 



in Figure |20| the effective energies at p = 1 and p = for j3 = 0.5 and A = 0.01. All 



results are tabulated in Tables 15 and [T(]. There is a general tendency for the vector 
boson mass to decrease, as the Higgs transition is approached, just as is observed for the 
scalar boson. The signal-to-noise ratio is too big to disentangle the dispersion relations 
in a way we did for the scalar ++ boson. Both one-particle dispersion relation and 
two-particle dispersion relation give consistent results within their error bars. Comparing 
with the masses extracted from the 'photon' operator Op at non-zero momenta (Tables 
|9| to 11) we find rough agreement. Therefore, the two operators - Oy and Op - seem to 



couple to the same state in the Higgs phase, the massive vector boson. This observation 
is in accordance with the perturbative analysis of the spectrum in the Higgs region, and 
confirms the early nonperturbative result fTof| obtained in the case of strong self-coupling 
(A = 3.0). 

In the Coulomb phase the situation is quite different from that in the Higgs phase. 
The correlators are very noisy and we do not have sufficient control of our signals to make, 
for example, an investigation of the existence of a bound-state like we attempted for the 
scalar boson. At non-zero momentum (there is no massless photon at zero momentum) 
the operator Oy should couple also to the photon. However, with rather large errors the 
effective energies for both p = and p = 1 are both quite large, ~ 2/a, and comparable. 
The effective energy for p = 1, in particular, is much larger than that of the photon, 
obtained from the correlations of Op. Evidently, the overlap with a photon state is very 
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Table 15: The effective energies for the 1 vector boson from unsmeared (normal) and 
smeared operators at time-slice 2 for A = 0.01. There was no signal for entries left blank, 
and none at all for (3 = —0.5. 



small. But due to the large errors on the effective energy, we can not make convincing 
arguments in favor of the existence of a vector bosonium in the Coulomb phase either. 
To settle the question of the excitation spectrum in the vector channel in the Coulomb 
phase would require much better statistics, or much better operators. 



6 Conclusions 

We have made a thorough study of the compact U(l) lattice gauge-Higgs model with com- 
pletely suppressed monopoles. We studied it both for comparatively small self-coupling 
(A = 0.01), as well as in the strong self-coupling limit (A = oo and A = 3.0). For strong 
self-coupling most computations were performed with A = 3.0, permitting us to compare 



our results with those of a previous study [||, 10 1 for the standard action. The computa- 
tions were performed for three different values of the gauge coupling (3: one at relatively 
weak gauge coupling, /3 = 2.5, one at stronger coupling, /3 = 0.5, and one at negative /3, 
corresponding to imaginary bare gauge coupling. 

We first mapped out the phase diagram of our model with monopoles suppressed. 
Contrary to the model with monopoles, it does not exhibit a confinement phase for strong 
gauge coupling. Nothing special seems to happen at [3 = 0, and both Higgs and Coulomb 
phases continue to negative [3 until about (3 = —0.7, where a first order transition into 
a frustrated phase seems to occur. The Higgs transition separating Coulomb and Higgs 
phase at positive (3 is clearly of first order for A = 0.01, similar to what happens in the 
model with the standard action M. Following Coleman and Weinberg |^], it could be 
interpreted as a radiatively induced spontaneous symmetry breaking. At strong self- 
coupling the Higgs transition becomes weaker and its order harder to determine. From 
our data we can not exclude that the transition becomes second order for sufficiently large 
A's, similar to the non-compact case E 



We measured T-dependent potentials from Wilson loops. In the Coulomb phase they 
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Table 16: The effective energies for the 1 vector boson from smeared operators at 
time-slice 2 for A = 3.0. There was no signal for entries left blank. 
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could be well fit with a T-dependent Coulomb potential (extracted from perturbative 
Wilson loops), even at the negative (3 probed, at least when the potential was extracted 
from smeared Wilson loops. The renormalized gauge coupling g\ grows with decreasing 
(3, and continues to rise rapidly at negative (3. The suppression of monopoles therefore 
does not seem to exclude a strong-coupling QED, but shifts it to the (unphysical) region 
9bare < 0- Since there the bare coupling is imaginary the transfer matrix is no longer 
positive and it is thus not clear whether the model makes sense physically. However, we 
would like to stress again that we detected no indication of a phase transition between 
the positive and negative j3 regions. 

In the Higgs phase we fitted the T-dependent potentials to T-dependent Yukawa 
potentials. This worked well for positive /3, and the extracted vector boson masses agreed 
reasonably with those obtained from conventional correlation functions. At negative f3 
the potential looked rather strange, and no good fits could be obtained. 

To study the spectrum of the model we calculated the correlation functions correspond- 
ing to the photon (at non-zero momentum), the scalar boson with quantum numbers ++ 
and the vector boson with 

Our results conformed to our expectations, in particular for the photon. In the 
Coulomb phase it was found to be massless (even at negative 0), and its energy fitted 
well the lattice dispersion relation. At k > k c (in the Higgs phase) the photon acquired a 
mass. The estimation of k c from the behavior of the photon is in a good agreement with 
that obtained from the Fredenhagen-Marcu order parameter. 

The measurements were more difficult for the ++ and 1 bosons, especially in the 
Coulomb phase. Both are reasonably well defined in the Higgs phase, as one would 
expect. Determination of the dispersion relation gave reasonably good evidence that the 
++ boson, the Higgs boson, is an elementary particle in the Higgs phase. In the Coulomb 
phase, however, the one-particle dispersion relation does not look very convincing, just 
like in the model with the standard action |To| . Within our large errors a two-particle 
dispersion relation, appropriate for a weakly bound state, did not work much better. 
Therefore, we do not think that the problem of the existence of bosonium, a bound-state 
between two scalar particles, in the Coulomb phase is settled. Much better statistics 
and/or better operators are needed to settle this question. 
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